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JH , The network of resonant bar detectors of gravitational waves resumed coordinated observations 



within the International Gravitational Event Collaboration (IGEC-2). Four detectors are taking 
part in this collaboration: ALLEGRO, AURIGA, EXPLORER and NAUTILUS. We present here 
the results of the search for gravitational wave bursts over 6 months during 2005, when IGEC-2 
was the only gravitational wave observatory in operation. The network data analysis implemented 
is based on a time coincidence search among AURIGA, EXPLORER and NAUTILUS, keeping 
the data from ALLEGRO for follow-up studies. With respect to the previous IGEC 1997-2000 
observations, the amplitude sensitivity of the detectors to bursts improved by a factor « 3 and 
the sensitivity bandwidths are wider, so that the data analysis wEis tuned considering a larger class 
of detectable waveforms. Thanks to the higher duty cycles of the single detectors, we decided to 
focus the analysis on three-fold observation, so to ensure the identification of any single candidate 
of gravitational waves (gw) with high statistical confidence. The achieved false detection rate is as 
low as 1 per century. No candidates were found. 

PACS numbers: 04.80.Nn, 95.30.Sf, 95.85.Sz 
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I. INTRODUCTION 

The search for transient gravitational waves (gws) re- 
quires the use of a network of detectors. In fact, the 
analysis of simultaneous data from more detectors at dif- 
ferent sites allows an efficient rejection of the spurious 
candidates, either caused by transient local disturbances 
or by the intrinsic noise of the detectors. Moreover, the 
false alarm probability of the network due to uncorre- 
lated noise sources at the different sites can be reliably 
estimated. 

The first long term search for bursts gw by a net- 
work of detectors has been performed by the five res- 
onant bars ALLEGRO, AURIGA, EXPLORER, NAU- 
TILUS and NIOBE, within the International Gravita- 
tional Event Collaboration (IGEC) [H. The search con- 
sisted in a time coincidence analysis over a 4-year pe- 
riod, from 1997 to 2000, and set an upper limit on the 
rate of impulsive gravitational waves as a function of the 
gw amplitude threshold of the data analysis 0]. How- 
ever, the overlap in observation time of the detectors was 
modest: three or more detectors were in simultaneous 
validated observation for 173 days, ~ 12% of the time, 
and two-fold observations covered an additional period 
of 534 days, ~ 36%. Moreover, since in the two-fold 
coincidence searches at the lowest amplitude thresholds 
some false alarms were expected, the IGEC 1997-2000 
observation was not able to discriminate a single gw can- 
didate from the accidental coincidences for most of the 
time. The target gw signals were short transients show- 
ing a flat Fourier component around 900 Hz, like pulses 
of 1 ms duration or oscillating signals with a few cycles 
of ~ 1 ms period. 

The same class of signals was targeted in the searches 
performed with the EXPLORER and NAUTILUS data 
in 2001 3 and 2003 0]. These searches, being based 
on two-fold coincidences, could not aim at the identifica- 
tion of single GW candidates. They addressed the study 
of a possible excess of coincidences taking advantage of 
sidereal time analysis. 

Subsequent searches for burst gws have been performed 
also by networks of intcrferometric detectors, which fea- 
ture a better sensitivity in a wider frequency bandwidth. 
In particular, the LIGO searches demonstrated a signif- 
icant improvement in sensitivity from the 2003 observa- 
tion to the 2005 observation [6|. Due to the shorter 
duration of these searches, an improvement on the limit 
set by IGEC at higher amplitudes on the rate of millisec- 
ond gw signals was not possible. However, In November 
2005, the LIGO observatory started its first long term 
scientific observation at its design sensitivity 0]. 

In 2004 four bar detectors resumed simultaneous oper- 
ation: ALLEGROji, AURIGA EXPLORER and 
NAUTILUS [nl, [ij A new long term gw search started 
under the IGEC-2 Collaboration, whose primary goal is 
to identify any single candidate of burst gw with high sta- 
tistical confidence. This coordinated observation is still 
running and targets to a broader signal class than the 



previous IGEC search, as for instance binary BH merg- 
ers and ring-downs [13| and longer transients recently 
predicted for Supernova core collapses [l3| . 

This paper is the first report on the IGEC-2 obser- 
vations and describes the results of the analysis of 6 
months of data, from May the 20*'* to November the 15*'' 
2005, when IGEC-2 was the only gravitational wave ob- 
servatory in operation. The AURIGA, EXPLORER and 
NAUTILUS data are actually used to search for gravita- 
tional wave candidates showing up as triple time coinci- 
dences. Due to a delay in the validation of the ALLE- 
GRO data, we agreed to use the data from this detector 
for follow-up investigations on possible signal candidates 
identified by the other resonant bars. 



II. CHARACTERISTICS AND GOALS OF THE 
IGEC-2 OBSERVATORY 



As for the previous IGEC search in 1997-2000, the de- 
tectors are aligned within a few degrees and so feature 
the same directional sensitivity at any time. The spectral 
sensitivities of the resonant bar detectors during 2005 is 
shown in Fig. [1] The minima of the noise power spectral 
densities are very close, within 1 2 x 10~^^l/^/Hz, as 
the four detectors share a similar design (i.e. cylindrical 
Al-5056 bar with a mass of ~ 2200 kg cooled at liquid 
He temperature, resonant transducers, similar mechan- 
ical quality factors, dc-SQUID signal amplifier). With 
respect to the previous IGEC-1 network 0], all detec- 
tors have been upgraded and exhibit now wider band- 
widths. EXPLORER and NAUTILUS were improved 
respectively in 2000 and 2002. Some modifications on the 
cryogenics apparatus and mechanical filters, new trans- 
ducers and dc-SQUIDs were adopted [l3|. The upgrade 
of AURIGA, completed in 2003, concerned most of the 
apparatus, from seismic isolation system [l^ to the read- 
out [l3l; in particular, a better coupling between the 
transducer and the signal amplifier was achieved by tun- 
ing the electric resonance of the signal transformer to the 
mechanical modes of the antenna and transducer and the 
signal amplifier is now based on a two stage dc-SQUID. 
The result was a very large increase in the detector band- 
width [13, [l^ . Additional upgrades of the room temper- 
ature suspensions during 2005 led to a significant im- 
provement of duty-cycle and data quality. ALLEGRO 
resumed operation in early 2004, after changing both the 
resonant transducer and the readout electronics Q. 

The primary scientific interest of the IGEC-2 observa- 
tions is the ability of identifying any single candidate of 
gravitational wave signal with high statistical confidence. 
Moreover, the results hereby presented refer to a period 
when only IGEC2 was surveying gws. 
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Frequency (Hz) 

FIG. 1: Typical strain noise spectral density (single-sided) 
curves of IGEC-2 detectors in 2005. From light gray to black: 
ALLEGRO, EXPLORER, NAUTILUS and AURIGA. All de- 
tectors share comparable minimum levels of noise spectra. 
The wider bandwidth of AURIGA includes the bandwidths 
of the other detectors. 



III. CHARACTERISTICS OF THE 
EXCHANGED DATA 

During the 180 days considered in this analysis, the 
AURIGA, EXPLORER and NAUTILUS detectors show 
a high duty cycle, see Tab. HI In particular the validated 
data of AURIGA, EXPLORER and NAUTILUS overlap 
by 130.71 days in three-fold coincidence (corresponding 
to 73%) and 45 days more are covered as two fold coin- 
cidences (about 25%). 

The noise stability of the detectors is remarkable, ei- 
ther if compared to past performances or to data analysis 
requirements. As shown in Fig. [21 the standard deviation 
of AURIGA noise shows a slow systematic dependence on 
the liquid He levels in the cryostats with peak-to-peak 
variations of the order of 10 %. Minimum noise levels 
of EXPLORER and NAUTILUS are higher by a factor 
of 2 in terms of equivalent amplitude of a millisecond 
gravitational wave burst. 

Each group validates its data and tunes its searches 
for gravitational wave candidates independently These 
analyses are based on linear filters matched to 5— like sig- 
nals. The algorithms implemented for the AURIGA filter 
and for the EXPLORER and NAUTILUS filter are dif- 
ferent and have been independently developed. In both 
pipelines, the filtered data stream is calibrated to give the 
reconstructed Fourier component H of the strain wave- 
form h[t) of a short {5 — like) gravitational wave burst 
at input. 

A candidate event is identified by detecting a lo- 
cal maximum in the absolute value of the filtered data 
stream: the occurrence time of the maximum and the 
corresponding amplitude are the estimates of the arrival 
time and of the Fourier component H of the gw h{t). 
Since these estimates refer to the 5— filter, they arc con- 
sistent only for gws of short duration with a flat Fourier 



transform over the detection bandwidth. For waveforms 
of colored spectral structure, the filter mismatch leads to 
non-optimal SNRs (Signal to Noise Ratios) for the can- 
didate events and to biases in their amplitude and time 
estimates. 

As an example, in the case of signals shaped as damped 
sinusoids with damping time r, the typical SNR re- 
constructed by 6 filters is > 80% of the SNR of the 
signal-matched filter for t < 10 ms, < 25 ms and 
< 50 ms for AURIGA, EXPLORER and NAUTILUS 
respectively. For such classes of waveforms, the imple- 
mented 5 — filters reconstruct the arrival times with rel- 
ative systematic errors < t/2. 

A cross validation has been performed on the differ- 
ent analysis pipelines. A sample day of raw data of EX- 
PLORER and of NAUTILUS was processed by AURIGA 
data analysis, using the same epoch vetoes, but with a 
different implementation of the data validation and con- 
ditioning. The comparison of the candidate events found 
by AURIGA and ROG pipelines show a good consistency 
for SNR > 5 with some unavoidable differences at lower 
SNR. 



TABLE I: Overview of the validated observation periods for 
the 180 days considered in this analysis. Off-diagonal terms 
show the two-fold coincidence times. 
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172.9 d 
151.8 d 
150.2 d 


158.0 d 
135.3 d 


155.0 d 
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FIG. 2: Noise vs time of AURIGA (black), EXPLORER (light 
gray) and NAUTILUS (dark gray) detectors. The ordinate 
shows la noise in terms of equivalent Fourier component H 
of the strain waveform h{t) of a millisecond gw pulse. 

Event lists of each detector are exchanged according 
to the protocol of the previous IGEC 1997-2000 obser- 
vations [igj . and include information on amplitude and 
time uncertainties and on the amplitude threshold used 
to select the events. 

A rigid time offset, chosen within ±10 s, is added to 
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TABLE II: Number of candidate events per each detector for 
some data selections considered in the network analysis. The 
leftmost column refers to the minimal thresholds, the central 
column to the event selection optimized for signals with com- 
parable SNR, the rightmost column to the selection optimized 
for signals with comparable H amplitudes, see Sec. IIVBI 



data cut AU 
EX 
NA 


SNR > 4.5 
SNR > 4.0 
SNR > 4.0 


SNR > 4.95 
SNR > 4.95 
SNR > 4.95 


SNR > 7.0 
SNR > 4.25 
SNR > 4.25 


AURIGA 


f869Il 


34598 


790 


EXPLORER 


489103 


29217 


245000 


NAUTILUS 


679775 


42028 


351375 



Amplitude H [ Hz ' ] 



FIG. 3; Amplitude distribution of the exchanged candidate 
events above the minimal thresholds: AURIGA (darker gray) 
SNR > 4.5, EXPLORER (lighter gray) SNR > 4.0 and 
NAUTILUS (gray) SNR > 4.0. The amplitude is given in 
terms of the Fourier component H of the h{t) waveform of a 
millisecond gw pulse. 



all the event lists prior to the exchange and is kept confi- 
dential, so that all the tuning of the analysis is performed 
without knowledge of the true coincidences. When the 
network analysis is completely defined, these confidential 
time shifts arc disclosed to draw the final results. This 
procedure, referred to as blind analysis, ensures an un- 
biased statistical interpretation of the results. 

The choice of the most suitable exchange threshold is 
left to each group. The thresholds used to select the 
exchanged events are SNR = 4.5 for AURIGA and 
SNR = 4.0 for EXPLORER and NAUTILUS. These 
arc considered as the minimal thresholds that allow the 
identification of a candidate event by each detector with 
reasonable confidence, according to the results of tests 
carried out with hardware and software injections. In 
particular, at lower SNR the timing uncertainty related 
to the candidate events increases rapidly. 

For the AURIGA events, the conservative estimates 
of the timing uncertainty (la) range from a maximum 
of 5 ms at the threshold to a minimum of ^ 0.5 ms at 
SNR > 10, as computed assuming S — like signals. For 
EXPLORER and NAUTILUS, the la timing uncertainty 
was conservatively set to 10 ms. 

The amplitude distribution of the exchanged events 
corresponding to the period of three-fold observations is 
shown in Fig. [S] The amplitude distribution of the AU- 
RIGA events is very close to that expected for Gaussian 
noise up to SNR ~ 5.5, while non Gaussian outliers are 
dominating at at higher SNRs. The number of candidate 
events above the minimal thresholds is listed in Tab. [TTl 
the mean event rate is ~ 45 /h for AURIGA, while it is 
larger for EXPLORER and NAUTILUS, - 129/h and 
~ 183/h respectively, due to the lower SNR thresholds. 



IV. NETWORK DATA ANALYSIS 

The network data analysis consists of a time coinci- 
dence search among the exchanged events. The coinci- 
dence time window is set accordingly to the same pro- 
cedure previously implemented within IGEC 1997-2000 
search [2| . Two events are defined in coincidence if their 
arrival times ti and tj are compatible within their vari- 
ances, erf and a'j : 

\U~t,\<k .Jaf + al (1) 

where k is set to 4.47, as in ref. 0. According to the 
Bycnaime-Tchcbychev inequality (see for example [20|). 
this choice limits the maximum false dismissal probability 
of the above coincidence condition to 5% regardless of 
the statistical distribution of arrival time uncertainties. 
For the three-fold coincidence search considered here, the 
same condition is required per each detector pair, leading 
to a maximum false dismissal probability < 1 — 0.95 ^ 
14%. The resulting coincidence windows are ~ 63 ms 
between EXPLORER and NAUTILUS and 45 - 50 ms 
when AURIGA is considered. In the previous searches 
performed between Explorer and Nautilus a fixed time 
window of 30 ms was adopted by the ROG group. This 
value ensured a low false dismissal in the case of delta- 
like signals considering the measured time response to 
excitations due to cosmic ray showers (ll| . 

Here we have neglected the effect of the propagation 
time of the gws among the different sites, since it is quite 
small, < 2.4 ms. Moreover, in case the signal duration 
is not small with respect to the coincidence window, the 
cited false dismissals are no more strictly ensured, be- 
cause the systematic uncertainties on the arrival time 
can be different in different detectors (see the previous 
section). 

The coincidence search is tuned to ensure a high sta- 
tistical confidence in case of detection of any single gravi- 
tational wave: 1 false alarm per century. To meet this re- 
quirement, we analyze only the 130.71 days of three-fold 
observation by AURIGA, EXPLORER and NAUTILUS, 
since the two-fold coincidence search cannot reach such 
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Number of coincidences per resampling 

FIG. 4: Histogram of the number of accidental coincidences 
per each resampling (black continuous line). Exchanged 
events with SNR > 4.0 for EXPLORER and NAUTILUS and 
SNR > 4.5 for AURIGA have been considered. 19'355'600 
independent off-source resampling of the experiment have 
been computed by shifting the time of two detectors within 
±11000 s in 5s steps, excluding the central region around 
the nominal zero-lag (see Sec. IIII|l . The histogram is well in 
agreement with the Poisson distribution with mean equal to 
2.16 (gray dotted line and shaded area), as = 11-3 with 12 
degrees of freedom corresponding to a p-value of 50.3%. 



a low rate of accidental coincidences without sacrificing 
too much on the sensitivity side. 



A. Accidental coincidence estimates 

The three-fold accidental coincidences has been inves- 
tigated with high statistics: about 20 millions indepen- 
dent, ojf -source, resampling of the counting experiment 
have been performed by applying relative shifts at the 
times of two detectors within ± 11000s in 5s steps. The 
resulting changes in the overlap time of the resamplings 
with respect to the actual observation time are negligible: 
the mean observation time of the resamplings is ~ 0.09% 
less than the actual observation time and the largest dif- 
ference is at most 0.4%. Fig. [4] shows the histogram of 
the counts of the accidental coincidences using the whole 
set of the exchanged events. The histogram is very well 
in agreement with a Poisson distribution of mean equal 
to 2.16 counts per observation time. 

Cross-checks on the accidental coincidences rate esti- 
mate have been performed with other independent al- 
gorithms and different choices of the relative shifts, giv- 
ing in all cases results well within the expected statis- 
tical fluctuations. A further check was pursued with a 
method based on an analytic estimate of the random co- 
incidence rates (see Appendix) . The results were in very 
good agreement with the values obtained with the time- 
shifts technique. 



B. Data selection 

In order to achieve the goal of 1 false alarm per century, 
a data selection is necessary to reduce by a factor w 600 
the accidental coincidences found on the exchanged data 
set. In general, the data selection has to be tuned with 
the aim of preserving the detection efficiency of the gw 
survey. In our case the balance between false alarms and 
detection efficiency has been addressed from first prin- 
ciples, since measurements of average efficiency during 
the observation time were not available. We decided to 
perform three searches based on different data selection 
procedures: 

• A) fixed and equal thresholds on the SNR of the 
events of each detector. Its motivation relies on 
setting a minimal comparable significance for the 
considered events as well as setting a similar events 
rate for each detector. Given the different spectral 
sensitivities, this search is more sensitive to col- 
ored signals that show the largest fraction of their 
power in the overlapping part of the bandwidths 
(e.g. 915 — 94:5Hz, see Fig. [1]) rather than in the 
remaining part of the AURIGA bandwidth. Such 
signals would produce similar SNRs in the S fil- 
tered data of all detectors. 

• B) fixed thresholds on the SNR of the events, but 
chosing different SNR thresholds for the detec- 
tors so that they correspond to comparable lev- 
els of absolute gw amplitude H in all detectors. 
This search is targeted to short signals which fea- 
ture a fiat Fourier transform within the AURIGA 
bandwidth and therefore appear at higher SNR in 
AURIGA with respect to EXPLORER and NAU- 
TILUS. It allows to use lower SNR thresholds for 
EXPLORER and NAUTILUS than the previous 
data selection procedure. 

• C) common absolute amplitude thresholds: same 
procedures used in the IGEG 1997-2000 search. 
The different data sets are selected according to 
a common gw amplitude Hi Q: the coincidence 
search is performed only during the periods when 
the exchange thresholds of all detectors were lower 
than Hi and considering only the events whose am- 
plitude is larger than Hi. This procedure is re- 
peated for a grid of selected Hi values. This search, 
as the previous one, is targeted to short bursts. 
Differently from the two previous procedures, this 
search not only selects the events, but also the effec- 
tive observation time as a function of Hi . Its main 
advantages are to keep under control the false dis- 
missal probability of the observatory and therefore 
to make possible an interpretation in terms of rate 
of gw candidates as well as a straightforward com- 
parison with the previous IGEG upper limit results. 

We decided to consider the union of these three 
searches, i.e. to perform one composite search made by 
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TABLE III: Number of accidental coincidences found on 
the 19355600 off-source resamplings per each data selection 
procedure (diagonal). The accidental coincidences found in 
common between different data selections are reported off- 
diagonal. Data selections B and C feature an evident correla- 
tion of their accidental noises. The false alarm rates of each 
trial taken separately are 0.396, 0.573 and 0.134 per century 
for A, B and C respectively. The resulting false alarm rate of 
the composite search, ^ U B U C is 1.01 per century. 



AU 


SNR > 4.95 


SNR > 7.0 


common 


EX 


SNR > 4.95 


SNR > 4.25 


search 


NA 


SNR > 4.95 


SNR > 4.25 


threshold 


data cut 


A 


B 


C 


A 


27368 






B 


515 


39507 




C 


147 


5177 


9280 



an "OR" of the three data selections procedures. This 
new approach simplifies the statistical analysis, since it 
takes care of the correlations expected in our multiple 
trials. In fact, any accidental coincidence occurring in 
more trials is counted only once and the expected over- 
all distribution of accidental coincidences is estimated by 
histogramming the union of the found accidental coinci- 
dences on off-source samples. 

Our tuning led to the following choices of data selec- 
tions: A) SNR > 4.95; B) AURIGA SNR > 7.0, EX- 
PLORER and NAUTILUS SNR > 4.25, C) common 
search thresholds 7?, = 1.2, 1.3, 1.4, 3.0 x lO-^i/iJz. 
Tab. |TT] reports the number of considered events of each 
detector for data selections A) and B). In particular in 
A) the event rate is similar in all detectors even though 
AURIGA features a wider bandwidth. Instead in B), 
the number of AURIGA events is a few hundred times 
smaller that that of EXPLORER and NAUTILUS. 

There is a significant correlation in false alarms be- 
tween data selections B) and C), which show a large frac- 
tion of common accidental coincidences. Tab. [HI] summa- 
rizes the numbers of accidental coincidences found for the 
three data selections on the same off-source resampling 
considered in Sec. IIV Al 

The final histogram of the accidental coincidences is 
plotted in Fig. O the probability of getting a non zero 
number of accidental coincidences is 0.00363 during the 
observation time, corresponding to 1.01 false alarms per 
century. The estimated la statistical uncertainty on this 
probability is 2 x 10"^. This uncertainty has been em- 
pirically determined by grouping the off-source samples 
in many disjoint subsets of equal size. The standard de- 
viation of the number of accidental coincidences in these 
subsets has been propagated to the mean, calculated on 
the entire off-source dataset. The resulting a is only 
slightly higher, by a factor ~ 1.4, than the one expected 
from a purely Poisson model. Independent checks with 
different pipelines and on different sets of off-source sam- 
ples limit the systematic uncertainty on the probability 
to < 1 X 10-"*. 
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FIG. 5: Histogram of the number of accidental coincidences 
found per each off-source resampling for the composite search 
made by the union of three data selection procedures (see 
Sec. IIV B]) . The histogram (black continuous line) is in agree- 
ment with a Poisson distribution with mean 0.00364 (gray 
dotted line and shaded area, = 0.06 with 1 degree of free- 
dom), which is taken as the reference distribution for the 
coincidences assuming that only accidental coincidences are 
present. The false detection probability is to 0.00363, corre- 
sponding to 1.01 false alarms per century. 



C. Plan of the statistical data analysis 

Before looking at the true coincidences in the on-source 
data set, we finalize a priori the plan for the statistical 
data analysis. Two steps are planned: the test of the null 
hypothesis and the setting of confidence intervals. 

We chose to reject the null if at least one triple coinci- 
dence is found in the on-source data set of the composite 
search. This corresponds to a significance of the test 
of 99.637% with a 3cr statistical uncertainty of ±0.006%. 
Therefore, if at least one coincidence is found, the collab- 
oration excludes it is an accidental coincidence with the 
above stated confidence: In fact, the rejection of the null 
points out a correlation in the observatory at the true 
time (i.e. not consistent with the measured accidental 
noise at different time lags). The source of correlation 
may be gws or disturbances affecting distant detectors 
(e.g. instrumental correlations). 

The final result on the estimated number of coinci- 
dences, related to any source of correlated noise or gws, 
is given by confidence intervals ensuring a minimum cov- 
erage, i.e. the probability that the true value is included 
in that interval. We decided to set confidence intervals 
according to the standard confidence belt construction 
of Feldman and Cousins [2l|. The noise model for the 
number of coincidences is the Poisson distribution shown 
in Fig. [5l To take into account its uncertainties, we con- 
sider the union of the confidence belts given by the mean 
noise ±3cr, i.e. 0.00364 ± 0.00006 events. Thanks to this 
low false alarm rate, the chosen confidence belt detaches 
from when at least one coincidence is found (provided 
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that the coverage is lower than the significance of the 
null hypothesis test). The final result cannot be easily 
interpreted in terms of gw source models, since IGEC-2 
is lacking a measurement of the detection efficiency. 

Any triple coincidence found would then be investi- 
gated a posteriori using also the data set of ALLEGRO, 
whenever possible. These follow-up results would be in- 
terpreted in terms of likelihood or subjective confidence 
by the collaboration and would not affect the significance 
of the rejection of the null hypthesis. The main goal of 
the follow-up investigation will be to discriminate among 
known possible sources, e.g. gravitational waves, elec- 
tromagnetic or seismic disturbances, etc.. An exchange 
of the raw data and gw transfer function would allow to 
implement more advanced network analyses, as searches 
based on cross-correlation. Additional complementary 
information could come from electromagnetic and neu- 
trino detectors as well as from environmental monitors. 



V. RESULTS 

Once the network analysis was tuned, the groups ex- 
changed the confidential time shifts necessary to recon- 
struct the on-source data set. This blind procedure 
makes the statistical interpretation of any result unam- 
biguous. 

No triple coincidences are found in the composite 
search described in the previous Sections and therefore 
the null hypothesis is not rejected. 

The upper limit set by the full search is given in terms 
of the number of detectable gravitational wave candi- 
dates, since the false dismissal of the composite search 
has not been directly measured for any model of gw 
source. According to the chosen confidence belt, the up- 
per limits are ~ 2.4 and 3.1 events at 90% and 95% cover- 
ages respectively. For a gw waveform with a flat Fourier 
transform over the bars bandwidths, the efficiency of this 
search is mainly contributed by the data selection B (sec 
sec. IIVB[) . In this case, according to back of the enve- 
lope calculations, IGEC-2 features a low false dismissal, 
< 0.1, at Fourier amplitudes > 2 x lO^^^ifz^^ for opti- 
mally oriented sources. 

Outside the planned composite search for gws, we 
checked a posteriori the number of coincidences among 
all exchanged events. Three coincidences were found, 
well in agreement with the expected Poisson distribu- 
tion of mean 2.16, as presented in subsection IIV Al All 
the events associated with these three-fold coincidences 
were at SNR close to the thresholds and therefore no 
follow-up investigation has been implemented for diag- 
nostic purposes. 



A. Comparison with IGEC previous results 

Using the subset of the current results relative to the 
data cut C, we can compute the upper limit on the rate of 
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FIG. 6: Comparison of the current upper limit with the pre- 
vious one set on 1997-2000 observations. The uninterpreted 
upper limits have been computed with the same methodology. 
The new result however uses only a subset (i.e. the C data 
selection, see Sec. IIlI[l of the composite search performed on 
2005 data. 



millisecond bursts as a function of the amplitude thresh- 
old. This upper limit is uninterpreted, i.e. it is set in 
terms of detectable gws, and is done mainly for com- 
parison with the previous IGEC 1997-2000 search Q, see 
Fig. El The new upper limit improves the old one at lower 
amplitudes thanks to the better sensitivity of current de- 
tectors. The current asymptotic rate, ~ 8Aevents/yr at 
95% coverage, is higher than in the previous search be- 
cause of the shorter observation time, but it is reached 
at much lower signal amplitudes. In fact, the current de- 
tectors feature much more stationary performances and 
the current search is free from false alarms, while the 
1997-2000 result was dominated at low amplitudes by 
two- fold observations, which typically show several false 
alarms per year. 

As a general remark, the main improvement of the 
current result is the capability of identification of any 
single candidate gw, while the previous upper limit was 
mostly contributed by coincidence searches with much 
higher false alarm rates. An additional improvement of 
the current search comes from the new data selections 
procedures (i.e. data selections A and B), which extend 
the target towards a broader class of signals. 



VI. FINAL REMARKS 

The IGEC-2 observatory is currently surveying for gw 
transients. Our results show that for a plain time coinci- 
dence search at least three-fold bar observations are nec- 
essary to identify any single candidate gw with satisfac- 
tory statistical confidence. The role of the resonant bar 
observatory is significative to search for signals occurring 



whenever the network of the more sensitive interferomet- 
ric detectors is not fully operative and therefore not able 
to issue an autonomous detection of a gw candidate. In 
fact, since the spectral sensitivity achieved by the LIGO 
instruments is better than a factor ^ 10 at the narrower 
bandwidths of the IGEC2 detectors, LIGO is nowadays 
able to perform surveys and, eventually, set upper limits 
at lower amplitudes and on a wider class of gws signals 
than IGEC-2 @. In this framework, IGEC-2 can col- 
laborate with the other observatories to extend the time 
coverage of current gw surveys and can contribute to the 
identification of rare gw events. In addition, if a candi- 
date gw will be identified by the interferometric obser- 
vatory, a joint investigation bar-interferometer could in- 
crease the information on the gw candidate, for instance 
on the signal direction and polarization amplitudes. To 
take the most from an hybrid bar-interferometer obser- 
vatory, the data analysis methodology should overcome 
the intrinsic limitations of a time coincidence search and 
exploit the phase information of the h{t) data streams 
provided by the different detectors, aiming at the solu- 
tion of the inverse problem for the wave tensor. Tests 
of such methodologies are ongoing using short periods of 
real data sets. 



APPENDIX A: ANALYTICAL ESTIMATE OF 
THE ACCIDENTAL COINCIDENCES 



The rate of accidental coincidences of the IGEC-2 
observatory has been empirically estimated by shift- 
ing the time of the detectors' data. The results have 
been checked by comparison with the following analytical 
method, based on the common assumption that the ex- 
changed events are Poisson point processes with a slowly 
variable rate. 



1. Analytical model 

The expected number of accidental coincidences Nacc 
in the simpler case of a constant time coincidence window 
±w and constant event rate is 



M-1 M 



(Al) 



i=i 



where M is the number of detectors, Tobs is the common 
observation time, and the number of events in the i*'* 
detector. 

In our case, the coincidence window w depends on the 
detector pair i,j and changes for the different AURIGA 
events (see eq.l and seclllip. Therefore, the above ex- 
pression for Nacc has to be modified as shown in the 
following. 

Given a time ti of an event of the detector #1, the 
probability that detector #2 has an event at a time ^2 



±(W,2- W,3) 




(W,2 + W,3) 



FIG. 7: P.d.f. of the variable x = t2 — is- The area of the part 
of the trapezium inside ±W23 (grey area) gives the probability 

of \x\ < W23- 



such that \t2 — ti \ < 11112 is 

Pl2 



2^N2 



(A2) 



and similarly for detector #3. These two probabilities 
arc independent, so that the probability that both occur 
is Pi2 ■ Pi3. When both occurs (necessary condition for 
a triple coincidence) we can find the distribution of the 
variable x = 1^ — 12 hy considering that the variables t2, 
have uniform distributions in the intervals ±wi2,±wi3 
respectively. Their probability density functions (p.d.f.) 
are then 



1 



2w 



(A3) 



with j=2 or 3, and their characteristic functions (Fourier 
transform of the p.d.f.) are 



$t,.(w) 



sin{uj Wij) 



OJ Wlj 

The characteristic function of the variable x 

<^^{U}) = $t3(w)$t3(w) 



(A4) 

t2 is 
(A5) 



and its p.d.f. F(x) is then given by the inverse Fourier 
transform of ^xii-o) 



Fix) 



1 

2^ 



+ 00 



(A6) 



yielding the trapezoidal shape shown in FiglT] and de- 
scribed by 



F{x) 



\W12 + Wl3 +x\ + \wi2 + Wi3 - x\ 
8W12W13 
\W12 - Wis +x\ + \W12 ~ Wi3 - 
8wi2Wl3 



(A7) 



where we have assumed Wi2 > W13. 

The Probability P23 that also detectors #2 and #3 are 
in coincidence, i.e. \t2 — is] < W23, is a fraction of the 
area of this trapezium: 



P23 — 



Pi: 



W23 
W12 ' 



for W23 < [wu - W13) 
for W23 > (wi2 + W13) 
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and in the intermediate range wi2 

Wi3 



Wl3 < W23 < Wi2 + 



which turns to eq lAll when all coincidence windows are 
equal to w. 



23 



2(«7i2Wl3 + W12W23 + W13W23) - (lUi2 + wf^ + 



23) 



4wi2Wl3 



(A8) 



The probability of a triple coincidence at each event of 
the detector #1 is given by the product P12 ■ P13 ■ P23- 
The number of accidental triple coincidences is obtained 
by further multiplying by Ni 



W12W13 
oet al.bs 



N1N2N3 



(A9) 



In our case the coincidence window is set from the tim- 
ing uncertainties cri.2,3 of the single events, according to 
eqUl Then, it is easy to verify that Wij is bounded be- 
tween the difference and the sum of the other two w's so 
that P23 is given by eq lASI The resulting analytical esti- 
mate for the accidental coincidences in case of constant 
event rates and different time windows is 



acc - T^NlN2N3{2{wi2Wi3 + W12W23 + W13W23) 



obs 



- W. 



23 



)} 



(AlO) 



r 



2. Implementation 

The common observation time of the three detectors 
has been divided in short sub-intervals with a dura- 
tion randomly chosen within a selected range, e.g. from 
~ 1 /2hour to ~ Ihour. The minimum and maximum du- 
ration must be chosen to meet the assumptions required 
by eq.AlO. In particular, the event rate should be sta- 
tionary, the coincidence window much smaller than the 



J 



average distance between events and the number of acci- 
dental (background) events much larger than the number 
of signal (foreground) events. 

Since the AURIGA events had a variable time uncer- 
tainty, we computed eq lAlOl for each of them using dif- 
ferent time windows. The prediction of Nacc is obtained 
by summing the result over all the AURIGA events in 
the J*'* sub-interval 



JVA„(i) 



NaccU) = ^NEAj)NNaij) ^ F (wiE , N) , w{A, , E) , w{Ak , N)) 



(All) 



k=l 



r 



where Tj is the interval duration, A''j 



Ex/Na/Au 



are the 



number of events of Ex, Na, Au, w{E, N) is the Ex-Na 
(fixed) coincidence window, w{Ak, E/N) is the coinci- 
dence window Au-(Ex or Na) computed with the at of 
the fc*'' Auriga event, and F{w,'w,w) is the combina- 
tion of windows in eq lAlOl The total result for the whole 



overlapping period is then obtained summing over all the 
sub- intervals. 

This procedure can be repeated with a different choice 
of the minimum and maximum intervals duration and/or 
with a different random initialization, in order to evaluate 
the fluctuations in the numerical value of the final result. 
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